Late universe dynamics with scale-independent linear couplings in the dark sector 
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We explore the dynamics of cosmological models with two coupled dark components with energy 
densities pA and pb and constant equation of state (EoS) parameters wa and wb- We assume 
that the coupling is of the form Q — H q(pA, Pb), so that the dynamics of the two components 
turns out to be scale independent, i.e. does not depend explicitly on the Hubble scalar H. With 
this assumption, we focus on the general linear coupling q — q + qA pA + qB pB, which may be 
seen as arising from any q{pA,pB) at late time and leads in general to an effective cosmological 
constant. In the second part of the paper we consider observational constraints on the form of the 
coupling from SN la data, assuming that one of the components is cold dark matter (CDM), i.e. 
wb = 0, while for the other the EoS parameter can either have a standard (wa > —1) or phantom 
(wa < — 1) value. We find that the constant part of the coupling function is unconstrained by 
SN la data and, among typical linear coupling functions, the one proportional to the dark energy 
density pA is preferred in the strong coupling regime, \qA\ > 1. Models with phantom wa favor a 
positive coupling function, increasing pa- In models with standard wa, not only a negative coupling 
function is allowed, transferring energy to CDM, but the uncoupled sub-case falls at the border of 
the likelihood. 

PACS numbers: 98.80.-k; 98.80.Jk; 95.35.+d; 95.36.+X 



I. INTRODUCTION 



The overall density of the observed universe, the 
growth of structures and their clustering properties can- 
not be explained by known forms of matter and energy 
[HQ- In addition, cosmic microwave background (CMB) 
anisotropy observations show that the total density is 
close to critical, so that the gap between known and un- 
known cannot be accounted for by curvature ■ Fi- 
nally, observations of type la Supernovae (SNe) [all 13 > 
baryon acoustic oscillations (BAO) [a [9(, and integrated 
Sachs- Wolfe (ISW) effect [H 113 tel1 us that the universe 
expansion is currently accelerating. To explain these 
facts cosmologists need to assume the existence of a dark 
sector in the theory [ijj], whose general properties have 
then to be tested against observations, i.e. with param- 
eters that have to be deduced from indirect evidences. 
One possibility is that the dark sector is accounted for, 
partly or in full, by a modified gravity theory 1 , while a 
more conventional approach is to assume that gravity is 
well described by general relativity, with the dark sector 
made up of an unusual energy momentum tensor. 

In the currently prevailing scenario, the dark sector 
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1 Sec e.g. 13] and other articles in the same special issue on dark 
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consists of two distinct contributions. One component, 
cold dark matter (CDM), accounts for about one third 
of the critical density [lj| and is needed to explain the 
growth of inhomogeneities that we observe up to very 
large scales, as well as a host of other cosmological ob- 
servations which goes from galactic scales, to clusters 
of galaxies, to redshift surveys. The other contribu- 
tion, dubbed dark energy, accounts for the remaining two 
thirds of the critical density, and is needed to explain the 
observed late time acceleration of the universe expansion 
[1, 0|- CDM can be modeled as a pressureless perfect 
fluid, representing unknown heavy particles, collisionless 
and cold, i.e. with negligible velocity dispersion. In its 
simplest form, dark energy consists of vacuum energy 
density, i.e. a cosmological constant A. Taken together, 
A and CDM make up for the the so-called concordance 
ACDM model [la Ha]- This simple model fits observa- 
tions reasonably well, but lacks a sound explanation in 
terms of fundamental physics, and a number of alterna- 
tives have been proposed. In general, dark energy can 
be modeled as a perfect fluid with an equation of state 
(EoS from now on) that violates the strong energy con- 
dition [17| . such that it can dominate at late times and 
have sufficiently negative pressure to account for the ob- 
served accelerated expansion. Scalar fields can also be 
formally represented as perfect fluids (see e.g. [l8[ and 
refs. therein). In a more exotic version, dubbed phan- 
tom energ y fl9l |20| , the EoS also violates the null energy 
condition [17l |. leading to the growth in time of the en- 
ergy density with the cosmic expansion. Finally, another 
rather radical alternative to ACDM is to assume a single 
unified dark matter (UDM), able to mimic the essential 
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features of ACDM that are necessary to build a viable 
cosmology. For example, in [21| we have considered ob- 
servational constraints on a UDM model with an "affine" 
EoS, i.e. such that the pressure satisfies the affine rela- 
tion P = P Q + ap with the energy density [22], [23[ . This 
model is a one parameter (a) generalization of ACDM, 
with the latter recovered for a = 0. There is no need to 
assume a-priori a A term in Einstein equations, because 
the EoS P = P + ap leads to an effective cosmological 
constant with fl A = -8wGP /[3H%(l+a)}. The problem 
is thus shifted from justifying a A term in Einstein equa- 
tions to that of justifying the assumed EoS: a possible 
justification of this affine model can be given in terms 
of scalar fields, either of quintessence or k-essence type 
[24] ]. This type of model escapes typical constrains on 
many UDM models [25[ (but cf. e.g. [26|]) because, for 
a given homogeneous isotropic background expansion, it 
allows multiple phenomenological choices for the speed 
of sound of the perturbations [27| • 

In models of the dark sector consisting of two compo- 
nents, dark matter and dark energy are usually assumed 
to interact only through gravity, but they might exhibit 
other interactions without violating observational con- 
straints [28[. Exploiting this degeneracy, here we depart 
from the standard scenario, and assume a cosmological 
model where the dark sector is made up of two coupled 
dark components, each described as a perfect fluid with 
its own constant EoS parameter w. This choice allows 
for the possibility that the observed evolution of the uni- 
verse, although reasonably well explained by the ACDM 
model, is actually due to the dynamics of two rather 
general coupled components, possibly alleviating the so- 
called "coincidence problem", fl\ as Qcdm, typical of 
the standard model [12 |. 

In this paper, our first aim is to characterise the dy- 
namics of our cosmological model with the two general 
coupled components, taking into account general forms 
of interaction, parameterized in terms of a late time func- 
tion Q linear in the energy densities, Eqs. (T5][ni). To this 
end we will use standard dynamical system techniques 
[2^ , l3(il ] , which are now rather common in the analysis of 
cosmological models, see e.g. [U HH, |33|, [13, IH, Hg| and 
[22I l23l 1371 [38j . To our knowledge, such an exhaustive 
analysis has not been carried out yet, although several 
sub-cases have been considered [H, HH H3, Ell, IH, 0, 0] . 
In our study, we restrict ourselves to the evolution of a 
homogeneous, isotropic cosmological background, leaving 
aside the question of what the effects of coupling could be 
in anisotropic models [23[ , or when general perturbations 
are present [45l Hfjj]. It is however worth noticing that, 
thanks to the particular form of coupling we choose, our 
analysis of the dynamics of the two components is valid 
in any theory of gravity, because is based only on the con- 
servation equations, and not on specific field equations. 

Secondly, as a way to gain some physical insight on 
the likelihood of some specific coupling models, we also 
explore the constraints on the predicted luminosity dis- 
tance modulus derived from type la Supernovae obser- 



vations, using a Monte Carlo Markov Chain (MCMC) 
approach. Needless to say, this is not intended as a full- 
fledged cosmological parameter estimation for these mod- 
els, but only as a first exploration of the parameter space 
to rule out those models which are manifestly in contrast 
with observations. This analysis requires the use of the 
Friedmann equation, hence general relativity is assumed 
as the valid theory of gravity. 

The paper is organized as follows. In Section II we 
study the general dynamics and solve the equations for 
the evolution of the coupled dark components; in Section 
III we focus on the cosmological effects of changing the 
parameters of the coupling term; in Section IV we derive 
constraints from the observation of type la SNe on some 
specific sub-classes of models; finally, in Section V, we 
present the main conclusions of our work. 



II. DYNAMICS OF DARK COMPONENTS 

A. Linear scale-free coupling 

In general relativity, assuming a a flat Robertson- 
Walker universe, the dynamics is subject to the Fried- 
mann constraint 
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where px is the total energy density of the various com- 
ponents. Beside baryons and radiation, pt includes any 
other component contributing to the dark sector, i.e. that 
part of the total energy-momentum tensor that in the 
context of general relativity is needed to explain the ob- 
served universe, in particular the CMB 0, [3], structure 
formation [lL |2| and the late time acceleration of the ex- 
pansion 0, H, S E3, El ■ The dynamics itself is de- 
scribed by the evolution of the Hubble expansion scalar 
H = a/ a, given by the Raychaudhuri equation 



H = -H A - 



(1 + 3wr)pi 



(2) 



This is coupled to the evolution equations for the energy 
density of each of the matter components contributing to 
Pt- Since H+H 2 = a/a, with a(t) the usual metric scale- 
factor (which we assume normalized to its present value), 
acceleration is achieved whenever wt = Pt / Pt < — 1/3, 
as it is well known. 

The standard ACDM model assumes two dark com- 
ponents: the pressureless cold dark matter (CDM), with 
wdm — 0, and the cosmological constant A with w\ = 
— I. CDM is needed to fill the gap between the baryon 
abundance and the amount of matter that is needed to 
explain the rotation curve of galaxies and structure for- 
mation in general, as well as to allow for a vanishing cur- 
vature model. In the context of general relativity, and 
under the Robertson- Walker homogeneus and isotropic 
assumption (see e.g. [47[ for alternatives), a cosmological 
constant A is the simplest possible form of dark energy 
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(DE) needed to generate the observed low redshift accel- 
eration. While this simple scenario is preferred from the 
point of view of model comparison and selection [2l[ , be- 
cause of the low number of parameters, from a theoretical 
perspective is oversimplified, and it is worth exploring al- 
ternatives, even if purely phenomenological. 

Here we shall consider two general coupled dark com- 
ponents with energy densities pa and p B . Since we want 
to introduce a rather general type of coupling, focusing 
our analysis on its effects, we shall assume the simplest 
possible form for the EoS of these two dark components, 
i.e. we will assume that the EoS parameters wa and w B 
are constant. On the other hand, we shall not a priori 
restrict our study to the sub-class of models where one of 
the two components represents CDM with, for instance, 
w B = 0. 

Due to the presence of the coupling, the two dark com- 
ponents satisfy the balance equations 

Pa + 3H(1 + w a )pa = Q (3) 

Pb + 3£T(1 +w B )p B = -Q. (4) 

Even assuming the linear form for the coupling Q given 
in ([5]) and ^ below, this model allows us to explore 
a large number of alternatives. Here we will focus on 
models for the homogeneous and isotropic background 
expansion, assuming that for those models that will fit 
current observational data it might always be possible 
to construct an appropriate perturbative scheme allow- 
ing for structure formation, for instance by assuming a 
vanishing effective speed of sound in one component. 

The coupled dark components pa and p B could be in 
principle be taken to represent DE only, i.e. they could be 
two extra dark components contributing to pr in |T]), in 
addition to CDM. Leaving aside this possibility, and ig- 
noring baryons and radiation as we will do in this Section, 
the sum of Eqs. ([3]|3]) gives the conservation equation for 
Pt = Pb + PA- A positive coupling term Q corresponds 
to a transfer of energy from p B to pa, and vice versa, but 
in general Q doesn't need to have a definite sign. 

An interaction term between two components has been 
considered several times in literature, starting from Wet- 
terich [H, [1| and Wands et al. Hz| in scalar field 
models, and has been analysed by Amendola in dark en- 
ergy model s (50l l5ll. l52l. 1531 . and for example recently in 

da, in so, EniiniinsiiEi. 

The coupling term Q can take any possible form Q — 
Q(H, pa, PB,t). Here we shall consider the case of an 
autonomous (t independent) coupling with a factorized 
H dependence 

Q=^Hq(p A ,PB). (5) 

As we shall see below, with this assumption the effects 
of the coupling on the dynamics of pa and p B become 
effectively independent from the evolution of the Hub- 
ble scale H . For this reason, we may call this a "scale- 
independent" coupling. Furthermore, with the decou- 
pling of the dynamics of the two dark components from 



that of H, the analysis of the next section is valid in any 
theory of gravity, because it is based on the conservation 
equations only: we don't need to use ([I])-©, i-e. the field 
equations of general relativity. Finally, we note that any 
coupling of this type can be approximated at late times 
by a linear expansion: 

q = qo + gap a + IbPb , (6) 

where <L4j<7b are dimensionless coupling constants, and 
go is a constant coupling term with dimensions of an en- 
ergy density 2 . In the following we shall analyse the dy- 
namics arising from this general linear scale-independent 
coupling. Obvious sub-cases are: q cc pr (qo — 0, qA — 
qs)', q oc pa {qo = 0, qs = 0); etc. We will come 
back to this in more detail in the next Section. Lin- 
ear couplings have been frequently analysed in literature 
(H El, \M, [50|,P, S HE HI) both for mathemati- 
cal simplicity, because they retain the linearity of sys- 
tem ([3]|4]) with no coupling, and because they can arise 
from string theory or Brans-Dicke-like Lagrangians after 
a conformal transformation of the metric. 



B. Analysis of the scale- free linear dynamics 

1. The linear dynamical system 

In order to proceed with the analysis of the dynamics 
of the dark components, let us change variables, using 
the total density pt = Pb + PA and the difference A = 
Pb — Pa- We also set 

w + = (w B + w A )/2 , w- = {w B - Wa)/2 , (7) 
q+ = (q B +q A )/2, q- = (q B - q A )/2 . (8) 

One reason for this choice is that ultimately the evolu- 
tion of pt is the one that governs the general expansion 
law through ([lj and (J2). In addition, thanks to the par- 
ticular form of the coupling ([5]) and assuming H > 0, the 
dynamics can be made explicitly scale- independent, elim- 
inating H by adopting N = In (a), the e- folding, as the 
independent variable. Then, denoting with a prime the 
derivative with respect to N, the system (011]) is trans- 
formed into 

p' T + 3pr(l + w+) + 3w_A = (9) 
A' + 3A{l + w + ) +3w-p T = -3(q+PT + q~ A + q ) (10) 

An effective EoS parameter u> e // is implicitly defined 
from Eq. ([9]): when u;_ = the two EoS coincide giving 



2 Strictly speaking, an expansion about today would lead to q = 
9() + Qa(pa — Pad) +<?s(ps — Pbo), but constants can always be 
rc-defined in order to put the coupling q in the form J6}. 
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rise to a constant w e // = w+ — wa and pr scales ac- 
cordingly, as a standard barotropic perfect fluid, but in 
general 



A 

W e f f = w + + w_ — 
Pt 



(11) 



changes with time. Notice that we can also define, using 
(JSJIHl) in ([3]), effective EoS parameters for the two compo- 
nents: 



WAeff = WA 
WBeff = Wb 



go + qBPB _ <1A 
2p A ' 2 

90 + QAPA <1B_ 

2p B 2 



(12) 
(13) 



From now on we will characterize the cosmological evolu- 
tion of any of the energy densities as standard/phantom 
behaviour. As mentioned in the introduction, stan- 
dard/phantom respectively correspond to an energy den- 
sity which is either a decreasing or an increasing function 
of time (the scale factor or the e- folding N) . The phan- 
tom behaviour arises in the presence of coupling from an 
effective EoS parameter < — 1, which corresponds to the 
violation of the null energy condition [T3| for that given 
energy density. Thus, it follows from (fTTj) and (|12j|13|) 
that we can have a phantom behaviour in the total en- 
ergy density p T as well as in one or both of the single 
components pa and ps, and that in principle the effec- 
tive EoS parameter of each of these can pass through the 
— 1 value, from phantom to standard or vice versa. On 
the other hand, we will also refer to constant parame- 
ters such as wa and % as having a standard/phantom 
value, respectively wa > — 1 or wa < —1, because the 
corresponding fluid would evolve in that way in the case 
of no coupling. 

We will also refer to an "affine" evolution. As said 
in the introduction, for an uncoupled component with 
energy density p this arises from an affine EoS of the 
form P = P Q + ap. Inserted in the energy conservation 
equation this leads to 



P = Pa + PQAia 



-3(l+a) 



(14) 



Therefore, starting from the Friedmann equations (Q][2|) 
with no cosmological constant term, the affine EoS and 
energy conservation lead to an effective cosmological con- 
stant pa plus an effective matter-like component with 
constant EoS parameter a (a barotropic perfect fluid) 
and today's density pqm (cf [2~il . I22I \2§L l2~ij for a de- 
tailed analysis of the cosmological dynamics arising in 
this case). As we will see, it turns out that there are 
solutions of the system (| MTU| that evolve according to 

©. 

In order to proceed with the analysis of Eqs. ([9l fT0|) 
using standard dynamical system techniques [29J, it is 
convenient to write it as 



where the phase-space state vector X and the constant 
C are 



Pt 
A 





-3<zo 



(16) 



and the matrix of coefficients J is given by 



— 3(1 + w+) — 3w_ 

-3(w_+q + ) -3(l + w+ + g_) 



(17) 



Fixed points, if they exist, are solutions X* of the equa- 
tion J X* + C = and, given that the system (fTS")) is 
linear, J is also the Jacobian of the system at these fixed 
points. These fixed points correspond to constant values 
of pt and A and in turn of pa and ps, that is to the 
emergence of an effective cosmological constant (when 
Pt =/= 0, see below). Every constant form of energy is 
indeed alike the cosmological constant A, and plays ex- 
actly the same cosmological role: when dominates the 
evolution of the background, it drives an exponentially 
accelerated expansion, with an effective EoS parameter 
close to —1. Therefore, in this Section we will focus on 
the analysis of these fixed points. 

Notice that - unlike the case with no coupling - there 
is no a priori guarantee from the equations above that 
PA and/or ps, as well as pt, will always be non-negative. 
However, one has to keep in mind that pt must be non- 
negative because of the Friedmann constraint (Q]). This 
means that if pt is vanishing for some value of N (a), 
then at that point the assumption H > 0, on the basis of 
which Eq. (fT5")) is derived, is violated, and the solutions 
of (|15p no longer correspond to solutions of the original 
coupled system of Eqs. © and (|3"ll4")) . 



2. Fixed points and stability analysis 

Even if system (|15p is linear, many different possibili- 
ties arise from the fact that it depends on five parame- 
ters 3 . There are two main cases, which we are now going 
to unfold. 

Case 1: det(J) ^ 0. In this case J -1 exists and there 
is a unique fixed point X* = — J _1 C. This can either 
be the origin in phase space, i.e. pt* — A* = 0, when 
go = (C = 0), or else this fixed point represents an 
effective cosmological constant, with pt* — Pa, A* = Aa 



X' = J X 



(15) 



3 Only four of them are independent: one can always rewrite Eqs. 
H9I10I I renormalising the parameters to any of them, but we don't 
want to assume that any of tu+, ui_ , , g_ , qo is non null. In 
particular, the value of go is irrelevant to the existence of the 
fixed points according to the eigenvalues of J, but physically its 
value is relevant, because it determines the effective cosmological 
constant 118> . 
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respectively given by 



Pa = 9 



w- qo 



Aa = -9 



det(J) ' 
(1 + w+) q 



dct(J) 



(18) 
(19) 



where 



det(J) = 9 [(l + w+)(l + W+ + q-) - W-(W- +q+)} , 

(20) 

with corresponding values for paa and pba- 

In order to analyse the stability properties of system 
(I15p at this fixed point, we now consider the eigenvalues 
of J. Given that det(J) ^ 0, there are two non-zero 
eigenvalues, given by 



A i 



tr(J) 



where 



tr(J) = -3 [2(1 +«;+)+<?-] , 

2 



D 



tr(J) 



- det(J) 



(21) 

(22) 
(23) 

(24) 



and the value of the discriminant D determines the three 
possible Jordan canonical forms of J. These three possi- 
ble cases, to be further clarified in the next subsection, 
correspond to D > 0, D = 0, D < 0, and are summarized 
below, with their sub-cases. 

Case la: real distinct eigenvalues (D > 0). There are 
three sub-cases: i) A_ < A + < 0, the fixed point is 
a stable node and both dark components have a stan- 
dard behaviour, with decreasing energy densities; ii) 
A+ > A_ > 0, the fixed point is an unstable node and 
both dark components have a phantom behaviour, with 
increasing energy densities; Hi) A+ > > A_, which re- 
quires det(J) < 0; the fixed point is a saddle and both 
dark components have a first standard phase with de- 
creasing energy density followed by a phantom phase. 
The connection between the phantom behaviour and sign 
of the eigenvalues will be rendered evident in the next 
Section, see Eq. (|2"5]) and related comments. 
Caselb: real equal eigenvalues (D = 0). In this case 
A + = A_ = Ao and the fixed point is an improper node, 
either stable (both dark components are standard) or un- 
stable (both dark components are phantom). In both 
cases the two dark components, as well as pt and A, fol- 
low a sort of affine evolution, with a modification term, 
see Eq. (j2"9"i 

Case lc: complex eigenvalues (D < 0). There are three 
sub-cases: i) if tr(J) = the fixed point is a centre; ii) 
if tr(J) > the fixed point is an unstable spiral; Hi) if 
tr(J) < the fixed point is a stable spiral. This last 
case is the most interesting, with px converging to an 



effective cosmological constant via a series of oscillations. 
Notice that this case has to be dealt with care, as in the 
past pr = (H = 0) at some point, and the time re- 
versal of system (1 WT0)) should be considered prior to that. 

Case 2: det(J) = 0. In this case if q ^ (C ^ 0) the 
system of linear equations J X» + C = is in general 
inconsistent and there are no fixed points. Alternatively, 
if qo = (C = 0) there is an infinite number of fixed 
points, each of them representing a possible asymptotic 
state depending on the initial conditions. These infinite 
number of fixed points is represented by a straight line in 
phase space, corresponding to a conserved quantity for 
system (TTSl) . It turns out that there are three possible 
combinations of the four parameters W- , q + , g_ that 
give det(J) = and they are shown in Table HI classified 
as Case 2a, 2b and 2c. For Case 2b, pt* = always, 
and either also A* = and then both pa* — Pb* — 
0, or A» ^ and then pa* = —Pb* = A*/2, so that 
one of the two is negative. Being the energy densities 
either null or negative, we can conclude that Case 2b 
corresponds to a non-physical situation. Case 2a and 
2c are more interesting, and are summarized in Table 
HIl In both cases, each fixed point on the line represents 
an effective cosmological constant. In particular, Case 
2a corresponds to wa = wb = —1, i-e. two cosmological 
constant-like components whose energy densities scale in 
a different way because of the coupling. In addition, it 
turns out that for this very peculiar case qo can be non- 
zero. 

Regarding the stability analysis, in Case 2 det(J) = 
implies that one of the eigenvalues is null, namely A_ = 
if tr(J) > (and vice versa) (cfr. Eqs. (HH) and ([2"3"|) V 
This implies that the total energy density, as well as the 
single dark components and A, follow the affine evolution 
(fT4"f ; we will comment further on this below Eq. (|2"5]) . In 
this case, the non zero eingenvalue A_/ + for each sub- 
case is given in Table |TJ Table [TT| gives instead the values 
the fixed points and the conditions for positive effective 
cosmological constants. 

Each of these fixed points is characterized by p\ cx 
Aa, which is equivalent to pba oc paa- Explicitly, for 
Case 2a the energy densities are related by pba = 
-pAAqA/qB, while for Case 2c p B A = -/9aa(1 + 
uu)/(l + wb)- Note that the same proportionality law 
holds for the fixed point of Case 1 (cf. Eq. (TT8lfT9]) ). This 
behaviour mimics that of scaling solutions [3l[|37l|. whose 
phase space typically admits fixed points where the con- 
tributions of the two fluids to the total energy density 
are constant (we will come back on this in Sec. IIIII) . It is 
easy to verify that at all these fixed points the effective 
EoS parameter (fTTj) has the value w e ff = — 1. 



C. An equivalent equation and its solutions 

While the classification of the various possible phase 
portraits of system ([MTU)) is best given by the represen- 
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Case 


parameters 


*+/- 


2a 


w- = , w+ - — 1 , Vg+ , Vg- 


-3fl- 


2b 


w- = 0, q- = -(1 + w+) , Vto+ , Vg+ 


-3(1 + w+) 


2c 


(l + w + )(l + w + +q^)-w± 

q+ = — , w_ , Vw+ , Vg- 


-3[2(l + iu + ) + g_] 



TABLE I: The three possible combinations of parameters giving det(J) = 0, cf. Eq. (|20|l . When det(J) = 0, one of the 
eigenvalues (121|) vanishes, giving a constant mode, and the other is given for each case in the third column here. Thus pT 
follows the affine evolutions (|14[1 . i.e. there is a constant mode and a power law (in the scale factor a) mode, like that of a 
barotropic fluid; A, pA and pB also have a a constant mode and the same power law mode. 



Case 


FPs 


Paa > 


Pba > 


Paa > , pba > 


2a 


PA = -^A A 


q+ > -q- 


q+ < q- 


A A > and -1 < q+/q- < , or A A < and < q+/q~ < 1 


2c 




R > -1 


R < 1 


A A > and -1 < 7? < , or A A < and < R < 1 



TABLE II: For Cases 2a and 2c, the fixed points and the conditions for positive effective cosmological constants, separated 
and combined, in both cases subject to the condition pa > 0. For Case 2a we have assumed go = 0, although this is not a 
necessary condition. For Case 2c we have defined R = (1 + w+)/w-. 



tation (|15p used in the previous Section, given that this 
is a linear system it is useful to consider the equivalent 
second order linear equation with constant coefficients, 
and its solutions. This helps interpreting the formalism 
of Sec. Ill Bl from a cosmological point of view. For this 
reason, we shall focus on the cases of main interest. 

First, notice that if ui- — then © decouples and 
implies 

p T = p c a^ 1+w +~> , (25) 

with p c an integration constant, so that the total energy 
density can either be standard or phantom. Eq. (|10[) can 
then be integrated, giving 

A = A a- 3 ( 1+w + +q -^ + A c a- 3 ( 1+W +^ + A, (26) 

where Ao is an integration constant, A c cx p c and A* cx 
go- This case is of limited interest, because in general 
the two components will have a negative energy density 
in the past or in the future. 

Assuming now u>_ 7^ 0, we obtain from (|9l ll0p a single 
equation for the total energy density 

p T - tr(J) p' T + det(J) p T = 9w-q , (27) 

with tr(J) and det(J) respectively given by Eqs. (j22|) and 
(I2TJ)) . It is noticeable that - when det(J) 7^ - this equa- 
tion admits a constant particular solution corresponding 
to the constant source term, given by the fixed point (fT8|) . 
Eq. |T5D with dct(J) 7^ is exactly this constant particu- 
lar solution, which is the effective cosmological constant 
term and it vanishes if = 0. 

When the eigenvalues (|2ip of the characteristic poly- 
nomial of Eq. (|2"T|) are distinct, the general solution can 
be written as 



When the eigenvalues are equal, A+ = A_ = Ao, we have 
the special Case lb of the previous Section and, denot- 
ing T = tr(J), we obtain: 

Pt = [pTi + PT2 ln(a)] a~ + p A . (29) 

Given that this evolution law is similar to the affine one, 
Eq. (H3}, except for the ln(a) correction term, and that 
the fixed point of this Case lb is an improper node, we 
can refer to this as an improper affine evolution. Here 
Pt+i Pt- and Ptu PT2 are integration constants related 
to the present values of the energy densities pa and pg 
of the two dark components; they are not independent 
in a flat universe, but related by imposing that the to- 
tal energy density has the critical value at present. The 
effective cosmological constant p\ is given by (fTS)) . In 
both cases (f2"5)l and (|29p the expansion of the universe 
is governed by this total energy density that enters the 
Friedmann equations (jTJS]) . The parameters (3± are re- 
lated to the eigenvalues A± by A± = — 3(1 + /3 T ). That is, 
f3 ± = f3 ± VD/3, with f3 = w+ + q_/2 = -tr(J)/6 - 1 
and D given by Eq. (f2U)) . The case (f2"9")) arises from 
D = 0, with A = -3(1 + A,) = tr(J)/2. 

It is interesting to note the role of the parameters 
(3± in (f28|) . assuming they are real. For a barotropic 
perfect fluid with constant EoS parameter w one has 
p cx a~ 3 ( 1+w \ Hence, (3± simply represent two constant 
effective EoS parameters that drive the evolution of the 
total energy density (|28| . In other words, the simple 
coupling (JSJIH]) of the two dark components with constant 
wa and in general produces a total energy density 
([28]) equivalent to that of two uncoupled fluids with con- 
stant EoS parameters /3±, plus the effective cosmological 
constant p\ 4 - 



p T = PT+a - 3(1+/3+) + PT -a-^ 1+ ^ + p A . (28) 



4 Mathematically, this is easily understood in terms of the dynami- 
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Apart from the asymptotic effective cosmological con- 
stant pa arising when q ^ and w_ ^ (Case 1 of the 
previous Section), it is also possible to have asymptotic 
effective cosmological constants when det(J) — (Case 
2 of the previous Section) if go = 0, corresponding to the 



vanishing of one of the eigenvalues (|21j) . i.e. A_/ + = 
and to /?+/_ = — 1. There are three cases, given in Table 
HI each with the corresponding non vanishing eigenvalue. 
In this cases Pt, A, pa and ps have the affine behav- 
ior (TT4")) . Hence the acceleration is guaranteed by the 
component of the total energy density for which an effec- 
tive EoS parameter /3+/_ assumes value —1. These fixed 
points seem less interesting however, because even when 
the non-zero eigenvalue is negative they are not general 
attractor of the dynamics (since the other eigenvalue is 
null): there is a different asymptotic non zero value of px 
for each possible initial condition. 

We have already commented above on Case lb, giv- 
ing (|29|) . We now consider the two other sub-cases of 
Case 1 of the previous Section from the point of view 
of the solution (f2"5)) . In both cases A, pa and pb evolve 
with the scale factor as pr in (j28|) , i.e. they are all lin- 
ear combinations of the two normal modes of the system, 
a -3(i+/3 + ) anc [ a -3(i+/3_) ) pj ug a cons t an t term. 

Case la. When both eigenvalues are real and distinct 
Pt scales as it would if we were considering two decou- 
pled components with EoS parameters (3+ and /3_, plus 
a cosmological constant. If det(J) < the fixed point 
is a saddle, i.e. unstable, with A+ > and correspond- 
ingly /3_ < — 1. Consequently pr scales as a mixture of a 
standard component and a phantom component. When 
dct(J) > either both /?+/_ > — 1, or vice versa, with 
corresponding standard/phantom behaviour. 

Case lc. If A± are complex conjugates, then f3± = /3 ± 
ii l /\D\/3 and the total energy density can be written as 



PT = a - 3{1+M {(p T - -p T +) sin y/\D\hi(a) 

+ {pt- +pt+) cos VTD|ln(a) }+pa (30) 



III. ANALYSIS OF SPECIFIC COUPLINGS 
A. Dynamics of density parameters 

Introducing an interaction between two fluids can lead 
to interesting solutions for the energy densities, like at- 
tractor points in the phase space where the contribu- 
tions of the two fluids to the total energy density are 
constants. In these points the value of the normalised 
energy densities depends only on the parameters of the 
model and, since they are attractors, they are reached 
from a wide range of initial conditions, thereby allevi- 
ating the coinciden ce p roblem. These are usually called 
"scaling solutions" [3ll [37j and are characterized by con- 
stant fractions of the energy density parameters, namely 
&a,b = Pa,b/{3H 2 ) (in units 8ttG = 1, c = 1). 

In order to analyse the dynamics of the system, let us 
define the new variables: 



PA 
3H 2 



y 



PB 

3H 2 



PA 

3H 2 



(31) 



where together with the coupled fluids we also include 
radiation to include the era when it's the dominating 
component, when initial conditions are usually set. Note 



that 



y 



£Ib and f2 7 are constrained by x 



y + fl-y = 1; z is the energy density parameter of the 
total effective cosmological constant, and we neglect the 
baryons contribution, which is always subdominant. The 
system (J3J4J) then becomes 



3(1 + w + - w-) + 2 



W_ 



(32) 



/ \ / \ det(J) 
(<1+ - Q-)x + (q+ + q-)y + — z 



9w- 



3(1 +w + +w-) +2 



(33) 



/ \ / \ det(J) 
(Q+ - Q-)x + (Q+ + Q-)V + — z 



9w_ 



„ H' 
H ' 



(34) 



containing an oscillating function that modulates the 
power law scaling. 

We will now investigate the dynamics of the density 
parameters. In Sec. lIVI we will present a MCMC analysis 
with SNe data. 



where 



cal system formalism of the previous Section: because the system 
1151 is linear, it can be globally (in phase space) written in Jor- 
dan normal form 29]. When the eigenvalues H21t X± are real 
and distinct this Jordan form is diagonal, i.e. the dynamical sys- 
tem separates into two normal modes, physically corresponding 
to two new effective uncoupled fluids with constant EoS param- 
eters f3±. 



w_ 
IT 



-1- -[x(l+3(w+ -«;_)) 

+y(l + 3(w++w-))+2(l-x-y)] (35) 



is a rewriting of the Raychaudhuri equation ^ for the 
Hubble expansion scalar. 

The fixed points, namely the points satisfying x' — 
y' = z' = 0, are presented in Table 11111 labeled by capital 
letters, together with the corresponding eigenvalues. To 
the best of our knowlegde, this is the first complete anal- 
ysis of the dynamics of a three components cosmological 
system where two of the barotropic fluids are coupled 
via a general linear coupling function of the form 



8 



The effective EoS parameters at each of the fixed points 
w e ff = Ptot/Ptot is also listed, where p tot = p A + p B + p 1 
and therefore w e ff = {w+ — w-)x + (w+ + w-)y + fi T /3. 

All the fixed points shown in Table HTTl exist for W- 7^ 0, 
when the EoS parameters of the two fluids are the differ- 
ent. As aforementioned, the only physically reasonable 
fixed point for system ([9l fl0|) corresponding to W- = is 
Case 2a, where det(J) = 0. 

The fixed points A corresponds to the radiation dom- 
inated era, while B, C and D represent epochs that are 
dominated by the two fluids. In particular, at the fixed 
point D the constant energy densities of x and y (paa and 
Pba) cause the accelerated expansion with w e ff = — 1. 
From the expression for x and y at this latter fixed point 
it is easy to see that y — —x(l + iua)/(1 + wb), which 
is exactly the proportionality that holds for Case 1 and 
Case 2c, as discussed at the end of Sec. Ill Bl this point 
is characterized by the final domination of an effective 
cosmological constant, either driven by qo (Case 1, see 
Fig. or not (Case 2c, see Fig. [3]). In the first case 
whenever \(3±\ < 1 D is always an attractor, while in the 
second case it is not because one of the eigenvalues is 
null. Notice that its existence is completely independent 
on q + and g_ . Whenever the system settles into the fixed 
points B or C the role of j3+ and j3- is exactly that of 
effective EoS parameters (see Table UlT|) which allows for 
phantom line crossing at late time (see Fig. [1]), i.e. line 
for which the effective total EoS parameter is w e f / = — 1 . 

In the following we will examine in more detail three 
special classes of the coupling function and in Sec. HVI we 
will make a first comparison of the models to the data 
using MCMC applied to type la SNe distance modulus. 



behaviour (|T4]) . or its improper modification ([29]) . is ob- 
tained in three cases: i) q = —2(wb ~ wa), which gives 
([2U)) : ii) q = — 2{1+Wb), that corresponds to /?+/_ = — 1 
(even for p\ = 0, i.e. go = 0, an effective cosmological 
constant is generated); Hi) wa = —1, where one of the 
two fluids is ab initio a constant term. Notice however 
that generally, because = + 2\/l5/3, models with 
(3 + = —1 and /?_ > — 1 are not feasible. In particular, 
the ACDM evolution is exactly recovered in case ii) for 
Wa = 0, that is if one of the fluids is dust; in case Hi) 
for qs = —2wb- 

The fixed point D is characterized by the domination 
of the constant part of the total energy density p\; along 
it, the values of x and y are both positive only if either 
wa or wb have phantom values, i.e. wa < — 1 or wb < 
— 1. This statement holds true also for models II and 
III. However, if wa < — 1 D is no longer an attractor, as 
A y = — 3(1 + w a) is greater than zero. On the other hand 
wb < —1 requires q > —2(1 + wg) to let the fixed point 
be an attractor: in this case q is positive. A strong and 
positive q corresponds to a transfer of energy from pa to 
the other fluid with wg < — 1. Therefore in order to fall 
at late time into the cosmological constant dominated era 
a fluid with a phantom EoS parameter wb must absorb 
energy from the other non-phantom fluid. It is worth 
stressing that the effective cosmological constant, i.e. qo, 
is somewhat redundant whenever the fixed point D is 
not an attractor (see Fig. [J) . In Fig. (fTJ) an example of 
this dynamics of the background is shown; the effective 
cosmological constant is not noticeable, since, after the 
evolution on the saddle point B, the system is trapped in 
the attractor point C. 



I. q+ 



II. q- = 



Imposing q + = g_ is equivalent to choosing qA = and 
qB = q', therefore among the range of possible couplings 
represented by (O we are restricting to the class of models 
where Q/H is proportional solely to the energy density 
of one fluid (in our case e.g. ps), and it reads 



H 



:(qpB + qo) 



(36) 



This assumption also includes models with q + = — g_ 
since the coefficients qA and qB can be cither negative 
of positive. In this case the dynamics is the same as for 
q + = q_, the roles of x and y being simply interchanged. 
We will refer to this subclass of models as model I. 

In this model y/~D is automatically real, since D = 
9(q/2 + wb — wa) 2 /4; as a consequence the scaling func- 
tion (|28[) always drives a power law expansion, with /3 + — 
q/2 + w B and /3_ = wa if flo > (i.e. (wb + wa + q/2) > 
0), vice versa if /3q < 0. Hence the total fluid ends up as if 
it was made up of: i) a component scaling as the original 
fluid pa with no coupling, ii) a, second component char- 
acterized by a new EoS parameter and Hi) an effective 
cosmological constant term p\. Moreover a pure affinc 



If =0 the resulting coupling function Q/H is lin- 
early dependent on the sum of the energy densities of the 
two fluids, approximately equivalent to the total energy 
density (these models have been examined for example 
in [40( and [58j) and is as follows 



Q 3 
H^2 {qPT 



qo) 



(37) 



With this assumption qA = qB = q+ = q and [3± = 
w + ± VD/3 where D = 9(wb — w A )(2q + WB — wa)/4. If 
q is positive these effective EoS are real for wb > wa or 
wb < W A — 2g, while if q is negative the same relations 
hold but with opposite inequality signs. We will label 
this model II. 

In this model the affine evolution is recovered for wb = 
(qwA+2wA + 2)/(q~2wA— 2), corresponding to/?- = — 1. 
In this case, which is indeed Case 2c of Section HH an 
effective cosmological constant arises even for p\ = 0. 
Again, because f3 + — /3_ + 2>/~D, models with + = — 1 
and (3- > — 1 are not feasible. From a cosmological point 
of view this means that a matter-like evolution cannot 
be generated together with a cosmological constant. The 
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Points 


X 


y 




w eff 




\ 

Ay 




A 











1 

3 


4 


1 - 3/3+ 


l-3/3_ 


B 


q--2w_+2VD/3 


(g_ -2to_ +2VD/3)(q + +2w_ +2VD/Z) 





P+ 


3(1 + Af) 


l-3C3 + -2 / 3_+8_)+ v /F + 


l-3(3+-25_+ 9 _)- v /F+ 


An]— 


4u>_(g++g_) 


2 


2 


C 


q_ -2w- -2\/Z>/3 


- 2u)_ -2v^D/3)(q + +2ii>_ -2-/D/3) 





P- 


3(1 + /?-) 


l-3(,fl_-2,3 + + q_) + ^F_ 


l-3(,fl_-2/3 + + 9 _)-^F_ 


4w_ 


4u)_(q + +q_) 


2 


2 


D 


1 + W— + ^+ 


1 — 10_ +li> + 


1 


-1 


-4 


-3(1+0+) 


-3(1 +/3_) 


2li>_ 


2^_ 



TABLE III: Fixed points of system (|32II34[) . the corresponding effective EoS and eigenvalues, where F± = 9qL /2 — 3g-(l ± 
V^S- 3w+) +9g+w_ +9w 2 _ -(l±2VS-3w+)(-l + 3w+). 
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FIG. 1: Upper panel: evolutions of the energy density pa- 
rameters Qa (thin solid line), Qb (dotted line) and fi 7 (thick 
solid line) for a model with q+ = q_ = 0.25; for compari- 
son, the dashed lines are the values of x and y at the fixed 
points B (thin short-dashed lines) and C (thick long-dashed 
lines). For this model the parameters are: £Iqa = = 0.5, 
WA = 0, wb = —1.5, /3+ = and f3- = —1.25. Lower panel: 
the total effective EoS parameter for the same model : Weff 
evolves from the value 1/3 in the radiation dominated era, ap- 
proaches the value in the matter dominated era and then 
asymptotically evolves toward a constant phantom value, in 
this case /3_ = —1.25. 



ACDM limit is achieved ifw A = (-1 + q ± ^q 2 + l)/2 
and wb = The evolution of the energy densities 

for a special choice of the parameters is illustrated in 
Fig. O the effective cosmological constant (fT5|) arises at 
late time, driving the acceleration, and p\ is caused by a 
non-zero go (^A ^ 0). 



1.5 
1.25 
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FIG. 2: Upper panel: evolutions for the energy density pa- 
rameters for a model with g_ = and q+ = —0.5; for com- 
parison, the dashed lines are the values of x and y at the fixed 
points B (thin short-dashed lines) and D (thick long-dashed 
lines). For this model the parameters are: QoA = = 0.5, 
wa = —1-1, wb = 0.2. Lower panel: effective EoS for the 
same model; for comparison, we plot the EoS parameter of 
the fixed point B, /3 + = —0.14. 

III. q+ = 

The subgroup of models with g+ = ( from now on 
model III) includes the couplings that are proportional 
to the difference of the energy densities A (for example 
recently analysed in [5f|). With this assumption g_ = 
qB = —qA = q and the discriminant D = 9(q 2 + (wb — 
uu) 2 )/4 is always positive, so that oscillating solutions 
([30]) are never permitted. The coupling function reads 

|=^A + go ). (38) 

As before, the affine expansion (jT¥J) may only be gener- 
ated if one of the two effective EoS parameters assumes 
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FIG. 3: Upper panel: evolutions for the energy density pa- 
rameters for a model with q+ — and q_ = —0.18; for com- 
parison, the dashed lines are the values of x and y at the 
fixed point B (thin short-dashed lines) and D (thick long- 
dashed lines). For this model the parameters are: Qoa = 0.5, 
Q.A = 0, u>a = —0.9, wb = 0. (dust). The EoS parameters at 
B are f3+ = —0.08 and 0- = —1. Lower panel: effective EoS 
for the same model. 

the value of the EoS of a cosmological constant, that is 
either 0+ = — 1 or /3_ = — 1. In particular if [3 + = — 1, 
/3_ = —1 — 2y/D/3 is always phantom. In this case 
none of the terms in Eq. (|28[) can play the role of mat- 
ter. On the other hand if /?_ = —1 (corresponding to 
K)4 = (-1 - q ± ^/l-q 2 )/2) 1 (3+ = -1 + 2y^D/3 is al- 
ways grater than —1, i.e. always standard. An example of 
this dynamics is shown in Fig. [3J where the effective cos- 
mological constant fTS]) arises at late time with no need 
of go, driving the acceleration (Case 2c). Then typically 
for wb = —q — wa — 1 we have that /3+ = and the 
ACDM model is recovered. 



IV. MARKOV CHAINS WITH SUPERNOVAE 

A. Methods 

Given the large number of parameters, the task of find- 
ing the minimum y 2 and mapping its distribution in the 
entire parameter space can be computationally expen- 
sive. To this end we adopt a MCMC. In this work we 
only want to test our models as a description of the ho- 
mogenous isotropic background expansion (regardless of 
perturbations), hence supernovas are ideal for this pur- 



pose. We use the 192 type la SNe distance modulus data 
set provided in [60j. In particular we want to see wether 
supernovae can qualitatively distinguish different kind of 
couplings, included what we called model I, II and III. 

Type la SNe light curves allow a determination of an 
extinction-corrected distance moduli, 

HQ=m- M = 5 log (d L /Mpc) + 25 (39) 

where d L = {L/AirF) 1 ' 2 = (1 + z) f* dz'/H(z') is the 
luminosity distance. We compare our theoretical predic- 
tions to the values of po with H 2 — 8ttG/3(pa + Pb + 
P- ( + Pb), where we account also for the baryon energy 
density pb- We fix the value of the dimensionless Hub- 
ble constant to be h = 0.72 [6l[ and the baryon energy 
density at present ttbh 2 = 0.02229 according to 0]- The 
smaller is the EoS parameter of a single fluid the later 
can be the domination era for this fluid. Hence, count- 
ing the role of (3± as effective EoS parameters, whenever 
f3± > a baryonic era might emerge at recent time. The 
absolute distance modulus M is intrinsically affected by 
uncertainty; therefore we treat it as a nuisance parameter 
and marginalize over it. 

The parameters that are representative of the 
models are {£loA,FlA,qA,qB,WA,WB}, or otherwise 
{Qoa, ^Aj q+i q~, w+t and, as functions of these, the 
two effective EoS introduced in Eq. (|28[) : f3+ and /3_. 
For the ensuing analysis it is worth reminding our clas- 
sification of models: I) model with a coupling function 
proportional to only one of the two energy densities; II) 
model with a coupling function proportional to the sum 
of the energy densities; III) model with a coupling func- 
tion proportional to the difference of the energy densities. 

We shall now focus our analysis on the case wg = 0, 
i.e. pb would represent standard CDM if it wasn't for the 
coupling with the DE component. 

B. Results: CDM - DE coupled models 

The first result we obtain is that Q\ is completely un- 
constrained, independently on which model we consider. 
This means that SNe are not sensitive to the constant 
term of the coupling. As we have seen in Sec. [Til the dy- 
namics of the system can easily generate the acceleration 
settling on fixed points D, where u> e // = — 1, even for 
Ha = (see Fig. [3]), or B and C, where the total energy 
density can also exhibit phantom evolutions. 

In Fig. 2] and O we present MCMC chains in a two- 
dimensional diagram ([qai <7b] on the right hand 
side). As said above, we consider a model where one of 
the two fluid represents a CDM component, i.e. wb = 0, 
a reasonable assumption considered all the other cosmo- 
logical probes pointing towards the existence of a form 
of cold dark matter (see e.g. 0), and we let wa assume 
three different values that characterize pa as a DE com- 
ponent (phantom-like behaviour is shown in the top pan- 
els, cosmological constant-like in the second row panels 
and non-phantom model in the bottom panels). 
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Note that the case where wa — and ps is DE can 
be easily derived from the previous one, corresponding in 
the diagram to a reflection with respect to the line qs = 
~qA- In fact, interchanging the two EoS and swapping 
the roles of the two energy densities, and applying the 
transformation (qA —> — qB, qB — > — qA, i-e. q+ —>■ —q+, 
g_ — » q~), one recovers the aforementioned model. 

In addition, the straight lines corresponding to models 
I, II and III are drawn, and diagrams of Fig. [5] are derived 
from the same choice of parameters as in Fig. |4] except 
for fl\ ^ (it is by eye easily verifiable that there is 
no dependence on Finally, the short-dashed curves 
represent the improper afhne evolution (|29p . while the 
short-dashed straight line represents affine models (fH|) 
with (3 +/ _ = -1 (Case 2c). 

As a first step we derived the unidimensional likeli- 
hood for {floA, q+, <?-}■ The best fit of the energy den- 
sity parameter for the three class of models presented in 
Fig. Q] and [5] is respectively f2 A = 0.63,0.65,0.76 with 
an error of 2 a = 0.1; this best fit does not change includ- 
ing f^A- In the diagrams £Ioa is therefore fixed to these 
best fit values. It is worth stressing that here we are 
not just analysing the typical models considered in liter- 
ature (namely I, II and III) but the results incorporate all 
the possible linear couplings, and, we might say, all the 
possible expansions at recent times of a generic coupling 
function Q (Eq. [5]) . Hence we are not interested in deriv- 
ing constraints on single parameters, a route that might 
be hard to follow with SN la in view the high number 
of parameters and their degeneracies. We instead want 
to see what kind of linear couplings are preferred by the 
data and provide a qualitative way to distinguish the type 
and the direction of the interaction. 

The first noticeable thing in the [<?a,9b] diagram is 
that the points lie almost on a horizontal branch of the 
diagram, close to the line representing model I, in partic- 
ular with Q oc pa- So if we allow the interaction term to 
be strong and move out of the weak coupling regime (i.e. 
|<7a,b| > 1), the most "frequent" linear coupling func- 
tion emerging from the chains is the one proportional to 
the DE density (pa)- In addition, strong couplings are 
favoured for positive value of qA (see Fig. 0] and O: the 
energy is transferred from dark matter to DE. Increas- 
ing the value of wa, that is moving from phantom- like 
values towards quintessence-like ones (going downwards 
in the right hand side column of fig. (01), this horizon- 
tal branch tends to negative values of qB- Models with 
a phantom wa show an increasing energy density with 
the scale factor, a, while for DE model characterized by 
wa > — 1 the energy density is diluted with the universe 
expansion: this second kind of models requires a lower 
transfer of energy from CDM to DE. Apart from a small 
spot in the origin of the axis (weak couplings), the cou- 
pling or type II does not seem to be favored by SN data, 
the effect increasing with higher values of wa, i-e. for 
non-phantom values. Another evidence that arises from 
diagrams Fig. [5] and [5] is that for non-phantom values 
of wa the uncoupled case (namely [qA,QB] = [0,0]) falls 



almost outside the border of the likelihood. 

Since today Pob — Poa we can say that the sign of the 
coupling function {Q ~ qAPoA + 9sPob) changes along 
the straight line qB = —qA (long dashed line): above 
this line the exchange term reverses the energy transfer 
from CDM to DE (i.e. positive Q), while below it is the 
opposite (negative Q). Again, the higher is wa the big- 
ger is the number of points that we can find below this 
line. Therefore for DE components with wa < —1 an 
exchange of energy from DE to CDM is less probable, in- 
dependently on the type of linear coupling. This reflects 
the fact that an increasing energy density (characteristic 
of phantom behaviour) favors more and more absorbing 
and positive DE couplings at present, while non-phantom 
values of wa seem to need a negative exchange term, most 
of all for weak couplings, to explain supernovae data. It 
is worth stressing that eventually the likelihood seems to 
exclude the uncoupled case. 

The connection between where the points lie in the di- 
agrams, i.e. the region favoured by the likelihood, and 
where the cosmological background evolution is affine 
is an interesting issue; this directly connects coupled 
DE models to an effective evolution of the total en- 
ergy density that is completely equivalent to a cosmo- 
logical constant plus a component with constant EoS 
parameter a, Eq. HU). If one looks at the left side 
diagrams of Fig. [5l a short-dashed curve and a short- 
dashed straight line are drawn on it. The former corre- 
sponds to the improper affine evolution (|29p , obtained for 
q- = ±^J —Aw-{q+ + w~). Hence the only affine models 
are those that correspond to the straight line for witch 
/?+/_ = —1 (for the model with wa — — 1 this coincides 
with the line representing model II, the only possibility to 
recover the affine evolution with no coupling) . For a DE 
model with a phantom wa the affine evolution coexisting 
with a non-zero fl\ is somewhat ruled out and, among 
the models indicated in the last Section, is more com- 
patible with a coupling function proportional to pa (DE, 
model I) and possibly to A (model III). For DE models 
with wa > — 1 the situation is different: the data seem 
to favor an affine evolution generated in models with a 
coupling function proportional to pb (matter, model I) 
and again model III. In addition, for DE models with 
standard wa an improper affine evolution together with 
a non-vanishing qo (J7a) is allowed, in a region where 
the coupling function shifts towards negative sign, thus 
representing a transfer of energy from DE to CDM. 



V. CONCLUSIONS 

We have analysed the dynamics of two coupled dark 
components represented by two barotropic perfect fluids 
characterized by constant EoS parameters wa and wb- 
We have assumed a flat, homogeneous and isotropic cos- 
mology and a general linear coupling between the two 
barotropic perfect fluids. This scale-independent cou- 
pling takes a linear form proportional to the single en- 
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FIG. 4: Coupling diagrams with two-dimensional likelihood for models with = 0. Apart from the short-dashed line 
that represents an affine evolution with /?+/_ = — 1, all the other lines are labeled with the corresponding type of coupling 
function (e.g. the solid line on the left side diagrams represents a coupling function Q oc pB (model I), while on the right side 
diagrams it represents Q oc pT (model II)). The energy density parameter at present is fixed at its best fit value, respectively 
n 0A = 0.63,0.65,0.76. 



ergy densities plus a constant term: any coupling of this 
type can approximate at late time a more general cou- 
pling function. We have studied the stability of the sys- 
tem and shown that an effective cosmological constant 
can arise both from the constant part go of the func- 
tion Q and from an effective cosmological constant-like 
EoS. We have also examined the dynamics of the energy 
density parameters, and evaluated the fixed points and 
the corresponding eigenvalues, for the most general form 
of linear coupling. We have then restricted the analysis 
to some specific linear couplings previously considered in 
the literature (model I, II, III). Since we are restricting 
to the background expansion and we have modeled the 



coupling function as a late time first order Taylor expan- 
sion, a comparison with distance modulus from SN la 
data appeared as our natural step further. We have pre- 
sented a MCMC analysis for a model with dark matter 
plus DE using the data set provided in [(5(| • Considering 
two representative specific values of the DE parameter 
wa, one standard (wa > —1) and the other phantom 
(wa < —1), we have condensed our results in coupling di- 
agrams, where the points arising from the MCMC chains 
are drawn together with lines for model I, II and III and 
for the improper affine (|29|) affine (fT4]) evolutions, the lat- 
ter including the ACDM model as a subcase. Couplings 
proportional to the DE density seem favored, mostly for 
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FIG. 5: Coupling diagrams with two-dimensional likelihood for models with Qa = 0.7. All the lines are labeled with the 
corresponding type of coupling function (e.g. the solid line on the left side diagrams represents a coupling function Q oc pb 
(model I), while on the right side diagrams it represents Q oc pr (model II)). The short-dashed line represents affine evolution 
with /?+/_ = —1 and the short-dashed curve represents affine evolution with /3+ = j3- . The energy density parameter at 
present is fixed at its best fit value, respectively Oqa = 0.63, 0.65, 0.76. 



strong couplings |g^| > 1. The total sign of the exchange 
term sets the direction of the interaction: models with 
phantom wa definitely prefer positive coupling, i.e. an en- 
ergy transfer from dark matter to DE. On the other hand, 
models with non-phantom wa not only allow for negative 
Q, but forces the uncoupled model to fall at the border 
of the likelihood. For further and stronger constraints 
more complementary data are required, like CMB spec- 
tra or matter power spectra. These observables necessi- 
tate an accurate relativistic perturbation analysis which 
is nor obvious neither uniquely defined in phenomenolog- 
ical coupled models as those considered here. Moreover, 



simplified observable that make no use of perturbation 
analysis, like the CMB shift parameter, can be strongly 
model-dependent and, although straightforward, should 
not be used in models where the evolution, even just 
that of the unperturbed background, detaches signifi- 
cantly from that of the ACDM model. These extended 
investigations can only be settled with future work. 
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